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Abstract 

A theory of dispersionless Andreev bound states on surfaces of time-reversal invariant uncon- 
ventional superconductors is presented. The generalized criterion for the dispersionless Andreev 
bound state is derived from the bulk-edge correspondence, and the chiral spin structure of the 
dispersionless Andreev bound states is argued from which the Andreev bound state is stabilized. 
Then we summarize the criterion in a form of index theorems. The index theorems are proved in 
a general framework to certify the bulk-edge correspondence. As concrete examples, we discuss 
(i) d X y-w&ve superconductor (ii) p^-wave superconductor, and (hi) noncentrosymmetric supercon- 
ductors. In the last example, we find a peculiar time-reversal invariant Majorana fermion. The 
time-reversal invariant Majorana fermion shows an unusual response to the Zeeman magnetic field, 
which can be used to identify it experimentally. 
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I. INTRODUCTION 



In unconventional superconductors, due to the sign change of the pair potential on the 
Fermi surface it is known that an Andreev bound state (ABS) via Andreev reflection 1-3 is 
generated at a surface of superconductors 4-9 . The ABS shows up as a zero bias conductance 
peak of tunneling spectroscopy 6-8 ' 10 , and through the study of tunneling spectroscopy of high 
T c cuprates 11-20 , it has been established that a dispersionless zero energy ABS 7 is generated 
for d^-wave superconductors. The existence of the ABS influences seriously interface and 
surface properties 21-26 . Great amount of unconventional quantum phenomena appear in 
quasiparticle tunneling 10 ' 27 ' 28 , Josephson effect 29-34 , spin transport phenomena 35 ' 36 , Meissner 
effect 37-40 and macroscopic quantum tunneling 41-43 . 

Dispersionless ABSs appear not only in spin-singlet d^-wave superconductors, but also 
in spin-triplet p^-wave ones 4 ' 5 ' 7 . However, it has been shown that the proximity effect of 
the ABS is completely different in these superconductors 44-48 . In the latter case, the ABS 
can penetrate into a diffusive normal metal (DN) attached to the superconductor, while it 
is not possible in the former one. The underlying physics is the existence of odd-frequency 
pairings 49 : Odd-frequency pairings are induced on boundaries of superconductors due to the 
breakdown of the translational invariance 50 . From the requirement of the Fermi statistics, 
the ABS in d^-wave superconductors is expressed by an odd-frequency spin-singlet odd- 
parity state, but the ABS in p^-wave superconductors an odd-frequency spin-triplet even- 
parity one 51,52 . Then the ABS in p^-wave superconductors can penetrate in DN since the 
odd-frequency spin-triplet s-wave component in p^-wave superconductors is robust against 
impurity scattering 40,49 ' 51 ' 53 ' 54 . 

Another important character of ABSs is that they realize Majorana fermions in condensed 
matter systems. For time-reversal breaking superconductors, ABSs with linear dispersion 
have been studied in the context of the chiral p-wave superconductor such as Sr 2 Ru04 55-58 . 
In this case, the ABSs are an analogous state to the edge modes of quantum Hall system 
(QHS) 59 ' 60 , which induce a spontaneous charge current along the edge 61 ' 62 . The main differ- 
ence between the edge state of QHS and the ABS is that the latter excitation is a Majorana 
fermion realizing a non-Abelian anyon, which can be used in a fault tolerant topological 
quantum computation 63 ' 64 . In addition to spin-triplet superconductors 60 ' 65-67 , now it has 
been known that spin-singlet superconducting states may support a non-Abelian anyon 68-74 , 
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and the various unusual transport properties have been explored 60 ' 65 ' 69-86 . For time-reversal 
invariant superconductors, the ABS with linear dispersion have been studied mainly in the 
context of non-centrosymmetric (NCS) superconductors 87-91 . One of the remarkable features 
of the NCS superconductors is that the superconducting pair potential becomes a mixture 
of spin-singlet even-parity and spin-triplet odd-parity ones 92 , and for CePtaSi 93 , the ABS 
and the relevant charge transport properties 66 ' 94-101 have been studied based on s + p-wave 
pairing. The resulting ABS has two linear dispersions analogous to helical edge modes 66 , in 
Quantum spin Hall systems (QSHS) 102-105 , thus instead of charge current, spin current is 
spontaneously generated along the edge The ABS forms a Kramers pair of Majorana modes 
called a helical Majorana edge mode, and several new features of spin and charge transport 
stemming from these helical Majorana edge modes have been predicted 66 ' 97-100 ' 106 . 

It has been found recently that Majorana fermions are possible also for ABSs with flat 
dispersion 107 ' 108 . Although s+p-wave superconducting states have been studied intensively, 
other types of NCS pairing symmetry may appear in strongly correlated systems. Micro- 
scopic calculations have shown that a parity mixed pairing state between the spin-singlet 
d x 2_ y 2-wave pairing and the spin-triplet /-wave one is realized in the Hubbard model near 
the half filling 109-112 . Moreover, a gap function which consists of the spin-singlet d^-wave 
component and spin-triplet p-wave one has been proposed, 113 as a possible candidate of 
superconductivity generated at heterointerface LaA10 3 /SrTi0 3 91 . Examining the ABSs in 
these NCS superconductors, 107 ' 108 we found a new type of ABS in the d xy + p-wave NCS 
superconductor 114 : Due to the Fermi surface splitting by the spin-orbit coupling, there 
appears a single branch of Majorana edge state with flat dispersion preserving the time re- 
versal symmetry. A similar ABS with flat band in NCS superconductors was also discussed 
in Ref. 115 . 

In this paper, we will address topological properties of the dispersionless ABSs on sur- 
faces of time-reversal invariant superconductors. For the ABSs with linear dispersion, it 
has been known that their presence is guaranteed by the topological invariance defined 
in bulk Hamiltonian. For instance, for full gapped two dimensional time-reversal break- 
ing (invariant) superconductors, a non-zero Chern number 116,117 (non-trivial Z 2 topological 
number 102 ' 103 ) ensures the existence of the Majorana fermions 59 ' 60 ' 118-124 , and for nodal su- 
perconductors, the parity of Chern number ensures the existence of Majorana fermions 74 . 
Furthermore, for spin-triplet superconductors, the ABSs on the boundary can be predicted 
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from the topology of the Fermi surface 67 ' 123 ' 125 . However, Majorana edge mode with flat 
dispersion and its relevance to the symmetry of Hamiltonian has not been clarified yet. It is 
an interesting issue to clarify the relevance of the topological property to the wave function 
of ABS with flat dispersion. 

The organization of this paper as follows. In section II, we will define a topological 
number starting from the bulk Bogoliubov de Gennes (BdG) Hamiltonian, and present a 
topological criterion for dispersionless ABSs. In section III, we discuss topological stability 
of the ABSs using the chiral symmetry of time-reversal invariant BdG Hamiltonians. Then, 
we propose index theorems for the dispersionless ABSs. In section IV, the relation between 
the sign change of gap function and topological criterion will be discussed. We will show that 
the sign change of the pair potential is directly relevant to the index theorem we proposed. 
As concrete examples, we will discuss the ABSs in spin-singlet d xy -wave superconductors 
and spin triplet p^-wave ones, respectively. By constructing the wave function of the ABS 
explicitly, the index theorems are confirmed. We also find that the topological structures 
of the ABSs are consistent with those of the odd-frequency pairings. In section V, we will 
apply our theory to the ABSs in noncentrosymmetric superconductors, in which a Majorana 
fermion with flat dispersion is realized. The wave function of the ABS for d xy + p-wave 
superconductors will be derived, and the relevant topological number will be discussed. 
It will be also found that the Majorana fermion has an anisotropic response to Zeeman 
magnetic field, which can be explained by a Z 2 topological number. In section VI, we will 
finally provide a general framework to certify the bulk-edge correspondence, and by using 
it the index theorems will be proved. In section VII, we will summarize our results, and 
discuss possible application of our theory to other systems. Throughout this paper, we use 
a convention of h — 1 except in Sec. VI. 

II. TOPOLOGICAL CRITERION FOR DISPERSIONLESS ANDREEV BOUND 
STATE 

We start with a general BdG Hamiltonian for superconducting states. 




(1) 
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with 

/ i(k) m , Affew 

\At(k) m , -£ T (-k) m ' 

where c\, a (cfc Q ) denotes the creation (annihilation) operator of electron with momentum 
k. The suffix a labels other degrees of freedom for electron such as spin, orbital degrees of 
freedom, sublattice indices, and so on. £{k) is the Hermitian matrix describing the normal 
dispersion of electron, and A(fc) the gap function which satisfies A T (— k) = —A(k). 

In the following, we suppose time-reversal invariant superconducting states. The time- 
reversal operation for electron is given by the antiunitary operator T in the form, 

T = UK, (3) 

where U is a constant unitary matrix acting on the suffix a of electron, and K the complex 
conjugate operator. K also reverses the sign of the momentum k. U satisfies UU* = —1 
since T 2 = — 1 for electron. Then, for time- reversal invariant superconductors, the BdG 
Hamiltonian %{k) satisfies, 

I U aa , \ 

0^(fc)0- 1 = n*(-k), = , (4) 

V U ^J 

which implies 

U£(k)U* = £*(-k), UA(k)U T = A*(-fc). (5) 
In addition, the BdG Hamiltonian H(k) has the following particle-hole symmetry, 



CU{k)C- x = -H*(-k), C 



( 



5 aa > 
5 aa i 



(6) 



Therefore, combining with (4) and (6), we obtain the so-called chiral symmetry 



! i // a- ! =o. i =-/ce - | ° lUL ' | . (?) 

-iU aa > 

which is the central ingredient of our theory. For later convenience, we perform here the 
following unitary transformation by which V is diagonalized 

UlTUr = | S °* ° V U T = ±( ^ ^ | . (8) 

-8 aa > ] v 2 \ -iU aa > -5 aa > 



Then the BdG Hamiltonian becomes off-diagonal, 

l4n(k)U r = \ ° q(k) I , (9) 

with 

q{k) =iS{k)W- A(fe). (fO) 

Now consider the ABS on a surface of the superconductor. The relevant topological 
number for the flat dispersion (or dispersionless) ABS is 

W(k\\) = --^ I dk^TU^ifydk/Hik)] 
Am J ' 



— Im 

2tt 



dk ± d k± In detg(fc) 



(11) 



where fey (k±) is the momentum parallel (perpendicular) to the surface we consider, and 
the integration is along the fcj_-direction with fixing the value of k\\ (Fig.l). For instance, 
for the surface perpendicular to the x direction, k± = k x and fey = (k y , k z ). 

For superfluids or continuum models of superconductors, the line integral in (11) is per- 
formed from k± = — oo to = oo. Far apart from the Fermi surfaces, we can regulate 
the gap function as A(fe) — > at = ±oo and neglect off-diagonal terms in £{k) without 
changing the physics. With this regularization, detg(fc) becomes identical at = ±oo. As a 
result, we can show that the topological number W(fc||) takes an integer value. On the other 
hand, for lattice models of superconductors, the line integral in (11) should be performed in 
the noncontractable closed loop C in the Brillouin zone as illustrated in Fig.2. In this case, 
the periodicity of the Hamiltonian with respect to a reciprocal vector G, H(k) = H(k + G), 
ensures the quantization of W(k\\). 

From the bulk-edge correspondence, we obtain the following simple criterion for the 
ABS 107 ' 108 ' 115,126 - 

• When the topological number W{k\\) takes a non-zero integer, a dispersionless ABS 
exists on the surface. 

We notice here that the resultant ABS has flat dispersion: This is because the topological 
number W(fc||) is nonzero in a finite region of k\\ since it cannot change unless the integration 
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(a) 



(b) 





FIG. 1: (color online) Surfaces of superconductors, (a) two dimensional surface, (b) 1 dimensional 
edge. The conjugate momenta of x± and xu are k± and kn, respectively. 



(a) 



(b) 













c 







FIG. 2: (color online) The integral path C in the first Brillouin zone, (a) two dimensional case, 
(b) three dimensional case. On the path C, the momentum k has a fixed fen. 



path intersects a gap node. From the above criterion, this implies that the zero energy ABS 
also exists in a finite region of fey. In other words, the ABS obtained from the above criterion 
has flat dispersion. On the other hand, if the integration path intersects a gap node, detq(k) 
becomes zero, thus a discontinuous change of W(fe||) becomes possible. (Note that at a gap 
node k , defH(k ) oc |detg(fc )| 2 = 0.) Therefore, we find that 

• the flat dispersion ABSs are terminated in a gap node in the surface momentum k\\ 
space. 

These results are confirmed in concrete examples in Sees. IV and V, and proved eventually 
in Sec. VI. 
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III. CHIRALITY AND TOPOLOGICAL STABILITY OF ANDREEV BOUND 
STATE 



In the previous section, we discussed the topological number constructed from the bulk 
BdG Hamiltonian H(k) which ensures the existence of the dispersionless ABS. In this sec- 
tion, we present an alternative argument on the stability of the flat dispersion ABS from 
a viewpoint of the boundary. This argument reveals a peculiar spin structure of the dis- 
persionless ABSs, which we call chirality, from which the scattering between them is much 
suppressed. 

In order to argue the boundary state, we perform the Fourier transformation of the 
BdG Hamiltonian V,(k) with respect to k±, and denote the resultant BdG Hamiltonian as 
H(x±, fcy), where x± is the conjugate coordinate of k±. Then consider the semi-infinite 
superconductor on x± > where the surface is located on x± = 0. The corresponding BdG 
equation is given by 

n(x ± , k\\)\u(x ± , fey)) = E(k\\)\u(x ± , fey)), (12) 
with the boundary condition 

|u(aj_L,fc||)) = 0, (13) 

at x± = 0. The zero energy ABS satisfies (12) and (13) with E(k\\) = 0. 
In our argument, it is convenient to use the following equation 

H 2 (x ± , k\\)\v(x ± , fey)) = E 2 (k l{ )\v(x ± , fc||)). (14) 

instead of the original BdG equation (12). As is shown in Appendix A, it is found that there 
exists one to one correspondence between the solutions of the BdG equation (12) and those 
of (14). In particular, the zero energy state \uq(x±, fey)) of (12) is exactly the same as the 
state \v (x±,k\\)} of (14) with E 2 (k\\) = 0. Therefore, we consider (14) instead of (12) in 
the following. 

The starting point of our argument is the chiral symmetry (7) of the BdG Hamiltonian. 
From (7), we have 

{H(x ± ,k\\),F} = 0, (15) 
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thus, H 2 (x±, feii) and T commute with each other, 

[H 2 (x±, fen), T] = 0. (16) 



Therefore, the solution \v (x±, fen)) of (14) can be an eigenstate of T at the same time. Then 

*\\) 



we denote the solution with the eigenvalue T = ±1 as \v^(x±, fen)). From the argument in 



Appendix A, we find the following properties of \v^ ± \x±, fey)). 

1. For a solution with nonzero E 2 (k\\) ^ 0, the state |i^ + ^(£Cj_, fey)) is always paired with 

~\x±, feii)) . In other words, for nonzero energy solutions, the number of the T = 1 
states is equal to that of the T = — 1 states. 

2. On the other hand, for zero energy solutions, the number of the T = 1 states is not 
always the same as that of the T = — 1 states. 

If we denote the number of the zero energy state with T = ±1 as N^, these properties 
imply that Nq + ^ — Nq ^ does not change its value against perturbation preserving the chiral 
symmetry. As illustrated in Fig. 3, some of the zero energy states might acquire non-zero 
energy by the perturbation, however, they always form a pair with opposite chirality T. So 
the difference Nq — Nq does not change as a result. This result implies that the existence 
of the zero energy state is robust against the perturbation once Nq —Nq~^ becomes nonzero. 

(a) 9 (b) 
E 2 E 







r=i r=-i 

o 



FIG. 3: (color online) Zero energy states of H 2 {x ± , fey), (a) A^ +) = 2 and Nq ) = 1. (b) = 1 

r(+) _jvrH 



and Nq ^ = 0. By perturbation, some of the zero modes in (a) might become non-zero modes as 



in (b), however, Nq ~Nq ' does not change. 



As was shown in the previous section, from the bulk-edge correspondence, we find that 
the nonzero W(fcn) implies the existence of the zero energy ABS. At the same time, in this 
section, we find that the nonzero 

N o +) ~ N t ] also 

ensures the robustness of the existence 
of the zero energy state. Therefore, it is naturally conjectured that these two quantities, 
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W(k\\) and Nq — Nq , should be equaled. Since there is a sign ambiguity to equal these 
two quantities, we have two possible equations in a form of index theorems, 

Wik^^N^-N^, (17) 

or 

W(k\\) = N^ +) - iV ( ->. (18) 

In the following examples, we will show that the conjecture of the index theorem (17) or 
(18) indeed holds. Then they will be proved eventually in Sec. VI. Here note that there are 
two possible choices of the surfaces of the superconductor which is perpendicular to x±, i.e. 
the surface of the semi-infinite superconductor on x± > or that on Xj_ < 0. It will be 
found that these two possible choices of the surface exactly correspond to the two possible 
equalities (17) and (18). 

IV. SIGN CHANGE OF GAP FUNCTION AND TOPOLOGICAL CRITERION 
IN SUPERCONDUCTOR PRESERVING S z 126 

We first consider the simplest case where the Cooper pair preserves spin in a certain 
direction, say S z . As is shown below, the gap function consists of a real single component, 
and the BdG Hamiltonian reduces to a 2 x 2 matrix in this case. Since the particle-hole 
symmetry (6) and the time-reversal invariance (4) have different forms in the 2x2 BdG 
Hamiltonian, it needs a special care to consider this case. In particular, the particle-hole 
symmetry (6) is not manifest in the 2x2 BdG Hamiltonian. Nevertheless, we find a chiral 
symmetry in the 2x2 BdG Hamiltonian in the following, and using it, we will present 
topological criteria similar to (17) and (18). 

In this case, it also has been known that a sign change of the gap function implies the 
existence of the zero energy ABS. 6 ' 7 Below, we will show that our topological criterion 
reproduces this result if the Fermi surface has a simple sphere-like shape. In addition, we 
consider a general Fermi surface, where these two criteria are not coincident with each other. 
It will be shown that our topological criterion excellently agrees with the details of the zero 
energy ABS, while the previous one is not. 

Let us consider a time-reversal invariant superconductor described by a single-band elec- 
tron. We also assume that the spin component of a certain direction is preserved. Without 
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loss of generality, we can select the preserved spin axis as z, and the gap function is given 
by 

{itp(k)a y , for spin-singlet superconductor 
id z (k)a z a y , for spin-triplet superconductor 

where ip(k) and d z {k) are real functions. Under this assumption, the Hamiltonian % reduces 
to 

^ = E( c !t> c ^)^x 2 (fc) I c t fct ] ' (20) 



with 



n 2x2 (k) = V } K ' . (21) 
A(Jfe) -e(Jfe) 



Here A(fe) is given by 



ip(k) for spin-singlet superconductor 
A(fc) = ^ . (22) 

d z (k) for spin-triplet superconductor 

Note that the BdG Hamiltonian is now reduced to the 2x2 matrix (21). 

As has been mentioned above, the particle-hole symmetry (6) is not manifest in (21). 
However, we find the reduced Hamiltonian (21) has the following chiral symmetry, 

{7,%x 2 (fc)} = 0, 7 = <7„. (23) 

Indeed, we can show that the chiral symmetry (23) is indeed a remnant of the original chiral 
symmetry (7). 

In a similar manner as Sec. Ill, we can define the topological number by using the chiral 
symmetry (23). By using the unitary transformation w 7 which diagonalizes 7 as 



%2x2(&) is recast into 



-1 / ' I i -1 / ' 



g*(fc) 
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with q(k) = —ie(k) — A(k). Now we introduce the topological number w(k y ) as 

W(ky) = -— J dk x tl [7^2"x2( fc )^^2x2(fc)] 



= — Im 

2tt 



dk x d kx In q(k) 



(26) 



Here the line integral (26) is performed in a manner similar to (11). As is shown in Appendix 
B 1, it is found that this integral can be rewritten as the following simple summation, 

= \J2 s ^Ak)} ■ 8gn[A(fc)], (27) 

e(fc)=0 

where the summation is taken for k x satisfying e(k) = with a fixed k y . The formula 
(27) makes it easy to evaluate the topological number w(k y ). Then from the bulk-edge 
correspondence, we can say that if w(k y ) ^ 0, there exists a dispersionless ABS on a surface 
of the superconductor which is perpendicular to the x-direction. 

Furthermore, by using an argument similar to that in Sec. Ill, it is found that the ABS 
with flat dispersion is an eigenstate of the chirality operator 7 again. Then, in a similar 
manner as (17), we can conjecture that the number of the zero energy ABS with chirality 
7 = ±1 satisfies 

w(k y )=n^-n^\ (28) 

or 

w(k y ) =n5 +) -nJ- ) . (29) 

These give a criterion for the ABS with flat dispersion. 

We notice here that our topological criterion of dispersionless ABS includes the criterion 
proposed previously. In the case where the topology of the Fermi surface is simple as 
illustrated in Fig. 4, it has been known that if the gap function satisfies 

A(k x ,k y )A(-k x ,k y ) < (30) 

then a zero energy ABS exists on the boundary perpendicular to x-direction. 6 ' 7 In other 
words, a sign change of the gap function with respect to k x — > —k x implies the existence of 
the zero energy ABS on a surface perpendicular to the x-direction. Our topological criterion 
reproduces this result correctly: The formula (27) leads to 

w(k y ) = \ [sgn[d kx e(-k° x , k y )} ■ sgn[A(-A;0, k y )} + sgn[d kx e(k° x , k y )} ■ sgn[A(A;0, ky)]] ,(31) 
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where (±k®,k y ) denote the intersection points between the integral path of w{k y ) and the 
Fermi surface. See Fig. 4. Noticing that sgn[dk x £(k®, k y )] = — sga[dk x s(— k®, k y )], we can 
rewrite this as 

w(k y ) = ±sgn[d kx e(-k° x , k y )} [sgn[A(-^, k y )] - sgn[A(^, k y )}] . (32) 

Thus the topological number w(k y ) becomes nonzero only when the gap function satisfies 
(30), which means that our topological criterion reproduces the previous one in this partic- 
ular simple case. 




FIG. 4: (color online) Fermi surfaces with simple topology in (a) quasi-one dimensional system, (b) 
quasi-two dimensional one, and (c) three dimensional one. The thick blue lines denote the integral 
path of w. For simplicity, we illustrate the integral path only in (a) and (b). For each case, the 
integral path gets across the Fermi surface only twice at k x = ±/c£. In (a) and (b), P and P' denote 
the intersection point (—k x ,k y ) and (k x ,k y ), respectively. 

We would like to emphasis here that our topological criterion does not merely reproduce 
the known criterion, but is more informative. First, the chirality of the zero energy ABS is 
determined in a manner consistent with the formulas (28) and (29). Second, our formula 
is also applicable to more complicated cases in which the previous criterion does not work. 
We will see them in Sees. IV A, IV B and IV D. 

A. d xy -wave superconductor 

Here we consider the d xy -weive superconductor where e(k) and A(fc) in (21) are given by 

e(*) = £-A*, A(*)=Ao^. (33) 
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V)(ky) 



Here A is a positive constant. From (27), the topological number w(k y ) is evaluated as 

1, for < k y < k F 
-1, for > k y > -kp , (34) 
0, for j/c^l > kp 

where fcp = y/2mfM is the Fermi momentum. Thus our topological criterion (28) implies the 
existence of the zero energy ABS for \k y \ < kp. 

Now, let us solve the BdG equation for the semi-infinite d xy superconductor on x > 
with the boundary condition \u(x = 0,k y )) = 0. For \k y \ < kp, the following zero energy 
ABS on x = is found 6 

\u (x)) = C\ 1 I e ik v y sm(k x x)e- x/ *, (35) 



— isgnk 



y 



where C is a normalization constant, k x = \Jk F — ky and = mA k y /k F . Since the ABS 
is an eigenstate of j(— a y ) with eigenvalue 7 = — 1 (7 = 1) for < k y < kp (0 > k y > —hp), 
it is found that = and = 1 for < k y < k F (n^ = 1 and = for 

> k y > —k F ). On the other hand, for \k y \ > kp, we do not find any zero energy ABS, thus 
= = 0. We summarize these results in Table I (a), which shows that the index 
theorem (28) holds in this case. 

If we choose the semi-infinite d xy superconductor on x < 0, the zero energy ABS on the 
surface at x = is given by 

\u {x)) = c[ 1 I e ik y y sm(k x x)e x ^. (36) 
\y isgnky J 

Thus and ^ are summarized as Table I (b). The index theorem (29) holds in this 
case. 



B. p^-wave superconductor 

Now consider the semi-infinite p^-wave superconductor on x > 0. The gap function is 
A(fe) = A k x /k with k = y/k^ and e(k) is the same as that in (33). For the p x -wave 
superconductor, we have 

1, for \k y \ < kp 

w(k y ) = <( , (37) 

0, for \ky\ > kp 
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(a) d xy -wave superconductor on x > 



ky 


n (+) 




n (+) 
' i o 




' i o 


w{k y ) 


< ky < k F 
> k y > —hp 
\k y \ > k F 



1 



1 




-1 
1 




i 

± 

-1 



(b) d xy -wave superconductor 


on a; 


< 


k y 


n (+) 




n (+) 




n ( - ] 


W{ky) 


< k y < k F 
> k y > — k F 
\k y \ > k F 


1 






1 



1 
-1 




1 

-1 





TABLE I: The number n[ ) ±) of the zero energy ABSs with the chirality 7 = ±1 for (a) the semi- 
infinite d X y-wme superconductor on x > and (b) that on x < 0. For comparison, we also show 
the topological number w{k y ) given in (34). The index theorem (28) and (29) hold in (a) and (b), 
respectively. 

from (27). Correspondingly, ii\k y \ < k F , we obtain the following zero energy ABS on x = 

\u(x,k v ))= ^ 1 ^e ik * y am(k x x)e- x/ S», (38) 
for the semi- infinite p^-wave superconductor on x > 0, and 

\u(x,ky))=\ 1 ] e ik y y sin{k x x)e x/ ^ , (39) 



for the semi-infinite p^-wave superconductor on x < 0. Here k x = \Jk F — ky and ^ 
mA /k F . It is also found that these solutions are the eigenstates of 7 with the eigenvalue 
7 = — 1 and 7 = 1, respectively. Thus and Uq ^ are summarized as Table II (a) and 
(b). We confirm the relations (28) and (29), respectively again. 
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(a) p^-wave superconductor on x > 



ky 


n {+) 




(+) (-) 
n o - n o 




\k y \ < kp 
\k y \ > k F 






1 



-1 



1 



(b) p x -wave superconductor on x < 


ky 


n (+) 




(+) (-) 

n ~ n 


W(*T/) 


\k y \ < kp 
\k y \ > kp 


1 







1 




1 





TABLE II: The number n ( ' of the zero energy ABSs with the chirality 7 = ±1 for (a) the semi- 
infinite p^-wave superconductor onx>0 and (b) that on x < 0, respectively. For comparison, we 
also show the topological number w(k y ) given in (37). The index theorem (28) and (29) hold in 
(a) and (b), respectively. 

C. Zero energy Andreev bound state and odd-frequency pairing 

In this subsection, we discuss the ABS in d xy -wave and p x -wave superconductors from a 
viewpoint of odd-frequency pairing. It will be shown that the odd-frequency pairing has a 
topological structure similar to that of the ABSs discussed in Sees. IV A and IV B. 

Odd-frequency pairing is the pair function (pairing amplitude) that changes a sign when 
exchanging the time coordinates of two electrons. 127 Thus the Fourier transformed pairing 
function is an odd function of frequency. Near the surface of superconductor, due to the 
breakdown of translational invariance, the pair potential acquires a spatial dependence lead- 
ing to the coupling between the even and odd-parity pairing states. From the Fermi-Dirac 
statistics, the induced pair amplitude at the interface should be odd in frequency. It has been 
established recently that zero energy dispersionless ABSs induce odd-frequency pairings at 
the surface of the superconductors. 49-52 

First, we consider a two-dimensional semi-infinite superconductor in x > where the 
surface is located at x — 0. As shown in the appendix D, the pair amplitude at the surface 
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of d xy -wave superconductor is given by 

/(fc.ig- 'g&jMMSJ (40) 

in the Matsubara representation. Thus the odd-parity odd-frequency pairing is realized at 
the surface. On the other hand, for p^-wave pair potential, we obtain 

= (41) 



Wn Vk 2 

which means the realization of even-parity odd-frequency pairing in this case. As well as the 
wave function of zero energy ABS of d xy and p^-wave superconductor in Sees. IV A and IV B, 
the factor sgnk y exists only for d xy -w&ve case. This factor decides the difference of the parity 
of induced Cooper pair. The difference of the parity results in a serious difference when we 
consider proximity effect into DN attached to superconductor 49 ' 51 ' 52 . In DN, only s-wave 
even parity pairing is possible. Thus, odd-parity odd-frequency pairing amplitude cannot 
penetrate into DN. This implies that ABS in d^-wave superconductor cannot enter into 
DN since it is expressed by odd-frequency spin-singlet odd-parity state. On the other hand, 
for p^-wave superconductor, ABS can enter into DN since it is expressed by odd-frequency 
spin-singlet even-parity state including s-wave channel 44-48 . 

Now we consider a two-dimensional semi-infinite superconductor in a; < 0. The corre- 
sponding pair amplitude at surface (x = 0) is given by 



isgn(k y )A \k x \\k y 
for spin-singlet d xy -wave superconductor and 



f(k X ,ky)=- ^^^*^^ (42) 



/(*xA) = -^W (43) 



lb \ l\r 

for spin-triplet j9 x -wave one, respectively. Comparing Eq. (40) [(41)] with (42) [(43)], we 
find the difference between them is the presence of — sign. The present — sign exactly 
corresponds to the different values of Hq and between case (a) and case (b) in Tables I 
and II. 



D. Dispersionless Andreev bound states in superconductors with multiple Fermi 
surfaces 

In the above two examples, we assume a simple Fermi surface obtained in (33). Here we 
consider superconducting states with more complicated Fermi surfaces. To realize multiple 
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Fermi surfaces, we consider a model in the square lattice. By taking into account the second 
next nearest-neighbor hopping, the normal dispersion e(k) is given by 

e(k) = — 2t[cos k x + cos k y ] — 2t'[cos 2k x + cos 2k y ] — //. (44) 

As illustrated in Fig. 5, if we choose the parameters t, t' and [i properly, multiple Fermi 
surfaces are realized in the first Brillouin zone. 



ky 







r 





FIG. 5: (color online) The multiple Fermi surfaces in the first Brillouin zone. We take t = if = 1 
and n = -2.5 in (44). 



For the gap function, we consider (a) d xy -wwe pairing with A(fe) = A sin k x sin k y , (b) 
Pa: -wave pairing with A(fe) = Aosin/ca; and (c) p x -wave pairing with A(fe) = Aosin2fc a; . 
Since the condition (30) is satisfied in all the cases in the above, the sign change criterion 
suggests the existence of zero energy ABS on a surface perpendicular to the x-direction. 
However, we find that zero energy ABSs do not always appear. 

To study the edge state, we consider the model above in the square lattice. The BdG 
Hamiltonian in the lattice space is given by 

% = "Hkin + H-s-i 

where i = (i x ,i y ) denotes a site on the square lattice, c] CT (ck )<r ) the creation (annihilation) 
operator of an electron with spin a at site i. The summation is taken for the 

nearest-neighbor (the second next nearest-neighbor) sites. H s is the pairing term depending 
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on the symmetry of the Cooper pair. For the d xy -wave pairing with A(k) = A sin k x sin k y , 
it is given by 

^ s = 4" ^2( C l+x+y,t C i,i + C i~x-y,t C l,l ~ C i+x~yAi,i ~ C i-x+y,t C i,l) + n - C " ( 46 ) 
i 

For the p^-wave pairing with A(fc) = A smk x , 

n ° = ^ YMA + ^ - c \A-*a) + h - c -> ( 47 ) 

i 

and for the p^-wave pairing with A(fc) = A sia2k x , 

H " = |f 5>lt c U,, - 4A-2x,i) + h.c (48) 

i 

Suppose that the system has two open boundary edges at % x = and % x = N x , and impose 
the periodic boundary condition in the ^/-direction. Solving numerically the energy spectrum 
as a function of the momentum k y , we examine the edge states. 

In Fig.6(a)-(c), we illustrate the energy spectra for each case. In spite that zero energy 
ABSs appear on the edges, the previous criterion does not explain the detailed structures. 
Indeed, in the cases (a) and (b), the dispersionless ABSs disappear near k y ~ 0, although 
the sign change condition (30) is still satisfied for both cases. Therefore, the sign change 
criterion does not work in these cases. On the other hand, the topological criterion does 
work in all the cases. From the formula (27), we obtain w(k y ) in Table III. Comparing w(k y ) 
in Table III with the zero energy states in Fig. 6, we find an excellent agreement between 
them. In particular, for (a) and (b), we find that w(k y ) = at k y ~ (i.e. ky > k y > —ky). 
Thus in contrast to the sign change criterion, our topological criterion explains the detailed 
structures of the dispersionless ABSs. 



V. NONCENTROSYMMETRIC SUPERCONDUCTOR 

In this section, we apply our topological criterion to time-reversal invariant superconduc- 
tors supporting multiple components of the gap function. For simplicity, we consider the 
superconductor in a single-band description. In a single-band description, the general BdG 
Hamiltonian for time-reversal invariant superconductor is given by 

, e(k)+g(k)-cT A(Jfe) 
H(k) = 1 \ yK ' K ' , (49) 

At(fc) -e(fe) + g(k) ■ a* 
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* 






Tj7 




FIG. 6: (color online) The energy spectra with open edges at i x = and = 100 (upper panels) and 
the corresponding gap functions in the first Brillouin zone (lower panels) for superconductors with 
multiple Fermi surfaces, (a) d xy -wave superconductor with A(fc) = Ao sin k x sin k y , (b) p x -wave 
superconductor with A(fc) = Aosin/c-j;, and (c) p^-wave superconductor with A(fc) = Aosin2A: x . 
We take t = t' = 1 and /j, = —2.5. Ao is chosen as Ao = 1 for (a) and Ao = 0.5 for (b) and 
(c), respectively. In the upper panels, k y is the momentum in the y-direction, and k y G [— 7r,7r]. 
The ABSs appear as zero energy states with flat dispersion. In the lower panels, the solid curves 
denotes the Fermi surfaces. The darker areas denotes the negative gap functions, the lighter areas 
the positive ones, and the gap functions vanish at the dashed lines. As a matter of convenience, 
we take the horizontal axis as the fey-direction in the lower panels. 



where e(fe) = e(— k) is the kinetic energy of electron measured from the Fermi energy, 
g(k) = —g(—k) the anti-symmetric spin-orbit interaction such as the Rashba spin-orbit 
coupling, A(k) the gap function A(fe) = iip(k)a y + id(k) ■ ao y . This Hamiltonian has the 
particle-hole symmetry (6) and the time-reversal invariance (4) with U = ia y . Thus it also 
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-1 
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- ^> /l^ ^> /by 











-fc^ > ky > — 7T 


-1 
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1 



TABLE III: Topological number w(k y ) for (a) the d xjr wave superconductor with A(fe) = 
Aosinfc^sin/cy, and the p^-wave superconductor with (b) A(fe) = Aosinfc x and (c) A(fe) = 
Aosin2/c x . We suppose that t = t' = 1, /it = —2.5 and Ao > 0. Here k l y (i = 1,2,3) are de- 
fined in Fig. 5. 

satisfies (7) with 

' °y \ / s 

(50) 



Oy 



For the dispersionless ABS for the semi-infinite superconductor on x > 0, the topological 
number (11) is given by 



W(k y ) = i-Im 



dk x dk x lndetg(fc) 



(51) 



with q(k) = [e(k) — iip{k)\a y + [g(k) — id(k)] ■ cra y . As is shown in Appendix B2, W(k y ) 
is evaluated as 



W(ky) = ~ J2 sgn[g(k)-d(k)-e(k)i;(k)]-sgn[d kx {e(k) 2 -g(k) 2 )}, (52) 

e (fc) 2 - S (fc)2=0 

where the summation is taken for k x satisfying e{k) 2 — g{k) 2 = with a fixed k y . 
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A. two-dimensional Rashba superconductor 



Now consider two-dimensional Rashba noncentrosymmetric superconductors. Here e(k) 
is given by (33) and g{k) = \{xk y — yk x ) where A is the coupling constant of Rashba spin- 
orbit interaction. Due to the Rashba spin-orbit interaction, the Fermi surface is split into 
two as illustrated in Fig. 7. The Fermi momenta for the smaller and larger Fermi surface are 
given by 

k x = -m\ + \/{m\) 2 + 2m//, (53) 

k 2 = m\ + \J (mA) 2 + 2m/i, (54) 

respectively. The spin-orbit interaction also mixes the parity of the gap function generally, 
so the spin-singlet component and the spin-triplet one coexist in the gap function 93,128,129 . 
The spin-triplet component d{k) is aligned with the polarization vector of the Rashba spin 
orbit coupling, d(k)\\g(k) 93 . Then, the triplet component is given by 

d{k) = A t f(k)(xk y - yk x )/k (55) 

with k = \fk? while the singlet component reads 

i>{k) = AJ(k) (56) 

with At > and A s > 0. The superconducting gaps are Ai = |Ai(fe)| and A2 = |A2(fe)| 
for the two spin split bands with Ai(fc) = (A t + A s )f(k) and A 2 (fe) = (A t - A s )f(k). 



k y 




FIG. 7: (color online) Fermi surfaces of two-dimensional Rashba noncentrosymmetric SC. 
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Using the formula (52), we evaluate the topological number W(k y ). We obtain 
W(k y ) = isgn[(A, - A s )f(-k 2 )\ - isgn[(A t - A s )f(k 2 )] 

-isgn[(A t + A s )/(-fei)] + isgn[(A t + A s )/(fc 1 )] (57) 

for \k y \ < ki, 

W(k y ) = isgn[(A, - A s )f(-k 2 )\ - isgn[(A, - A s )/(fc 2 )] (58) 

for k\ < \k y \ < k 2 , and W(k y ) = for k 2 < \k y \. Here, ±fcj = {±<\Jkf — k 2 , k y ) (i = 1, 2). 

From Eqs.(57) and (58), it is found that W(k y ) is nonzero only when f(k) is an odd 
function with respect to k x . Therefore, for example, for an s + p-wave or a d x 2_ y 2 + /-wave 
Rashba superconductor, where f(k) is given by f(k) = 1 or f(k) = (kl—ky)/k 2 , respectively, 
W{k y ) becomes zero. Consistently, we find that they do not support the dispersionless ABSs. 
On the other hand, for a d xy +p-wave superconductor, we have a dispersionless ABS. 



1. d xy + p-wave superconductor 

For the d xy + p-wave superconductor, we have f(k) = k x k y /k 2 . Thus, W(k y ) can be 
nonzero value. From Eqs. (57) and (58), we obtain 



W(k y ) = { 



for A s > A t , and 



W(k y ) = { 



2sgnk y , for \k y \ < k-y 
sgnk y , for k\ < \k y \ < k 2 , 
0, for \k y \ > k 2 

0, for \k y \ < ki 
—sgnky, for k\ < \k y \ < k 2 , 
0, for \k y \ > k 2 



(59) 



(60) 



for A t > A s . 



To confirm the relation (17), let us consider the two-dimensional semi-infinite d xy +p-wave 
Rashba superconductor oni>0 where the surface is located at x — 0. The corresponding 
wave function is given by 98,107 

\u(x, k y )} = [cfipt ex.p(iqf x x) + c^i" ex P(-^i^) + 4^t exp(ig^x) + cJ^J exp(-iq^ x x)] 
xexp(ik y y), (61) 
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± _ k± k m E 2 -[A m (kf {2))] 



with K(2) X = K(2) X = \f k l(2)- k y for \K\ < k H2) and fc+ (2)x = -A; 1(2)x = i^/fcg - fc 1(2) for 
|^| > fci( 2 ), and fe^ 2) = (±^( 2 ) x ,^j/)- Here = 1 ' 2 ) are S iven b y 



with 



T ^ = (l, -ia x ±, ia^Tii, T 1± ) , (63) 

T i>t = (i, ««2± ? *«2± r 2±, -r 2 ±) (64) 

Al(2)(fcf (2) ) 

1 1(2)± - , _ l° 5 J 

^±^ 2 -[A 1(2) (fef (2) )P 

and ai(2)± = (±k^ 2 ^ x — ik y )/ki( 2 y E is the quasiparticle energy measured from the Fermi 
energy. For E 2 < [Ai( 2 )(fc^ 2 J] 2 , the branch of the square root in Eqs.(62) and (65) is chosen 
so as the wave function (61) is normalized (i.e. \u(x, k y )) — > as x — > 00). When E = 0, we 



find 



-isgnky, for |A^| < fci 
ri± = S , (66) 

±isgak y , for > fci 



and 

-isgn(A t - A s )sgak y , for \k y \ < k 2 
±isgn(A t - A 8 )sgnhy, for \k y \ > k 2 

Thus, it is found that ipi± and ip 2 ± are eigenstates of T in (50), 

( -agnhyipf, for < fci 
isgnA^^ , for Ifcyl > fci 



(67) 



(68) 



. sgn(A t - A s )sgnfcj,^ 2 , for |fc y | < fc 2 

r W = s + • ( 69 ) 

=Fsgn(A t - A s )sgnkyil>2 , for \k y \ > k 2 

In Table IV, we summarize the chirality of ipf^y 

To construct the dispersionless ABS, we put the boundary condition (13) on the wave 
function above. Then, we have 

C +Vi + + cr^r + + ^^ 2 ~ = 0. (70) 
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The dispersionless ABS is obtained if there exist non-zero cfs satisfying (70). Applying the 
chiral projection operator P± — (1 ± T)/2 on the both side, we can divide (70) into the 
sector with T = 1 and that with T = —1. For example, when A s > A t and k 2 > k y > ki, we 
obtain 



in the T = — 1 sector. Then, we solve (70) in each chiral sector. In the above example, we find 
that in (71) is identically zero, but there exists a single non-zero solution of (c^,c^,c^) 
satisfying (72). To see this, let us notice that only two components of ipf are independent 
when tpf is an eigenstate of T. Hence, we obtain two independent linear equations from 
(72). Solving these linear equations, we have a unique solution (ci,c 2 ,c 2 ) of (72) up to a 
overall normalization factor. This result means that no ABS exists in the T = 1 sector but 
a single ABS in the r = — 1 sector. Thus, Nq + ^ = and Nq ^ — 1. In a similar manner, 
we can solve (70) for other k y s\ In general, we find that if one of the chiral sectors consists 
of three (four) wave functions, we have a single non-trivial solution (a pair of non-trivial 
solutions) of (70). In Table IV, we summarize the number of dispersionless ABS in 
each sector obtained in this manner. Comparing with the topological number W(k y ) given 
in (59) and (60), we find that our result agrees with the index theorem (17) exactly. 

In a similar manner, we can also construct the zero energy ABSs on the surface of the 
two-dimensional semi-infinite d xy +p wave Rashba superconductor on x < 0. In comparison 
with the ABS of the semi-infinite superconductor on x > 0, the chirality of the ABS is found 
to be opposite. Thus we find that the index theorem (18) holds in this case. 

We notice here that the zero energy bound state for k 2 > \k y \ > k\ is a Majorana 
edge mode. The wave function for the zero energy edge state ty m (k y ) can be written as 
T ^m{k y ) = (u 1 (k y ),U2(k y ),v 1 (ky),v 2 (ky)) where 



cfi/jf = 



(71) 



in the T = 1 sector, and 



(72) 



Ul (k y ) 



u 2 (k y ) = iavi{ky) 



icrv 2 (k y ) 



(afi - gi/zj exp(ik y y - if) 
_ (fi + Mi) exp(ik y y - if) 



(73) 



(74) 
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with a = (k y — ^Jky — kf)/ki, j3\ = (ak y /k 2 + 1), $2 = (a + k y /k2) and a = sgn(k y ). The 
functions /i and f 2 decays exponentially as a function of x and are even function of k y . The 
Bogoliubov quasiparticle creation operator for this state is constructed in the usual way as 
^(k y ) = u 1 (ky)c\(k y )+U2(ky)cl(k y )+v 1 (k y )c t (-ky)+v 2 (ky)ci(-k y ). Since ui(k y ) = v{(-k y ) 
and U2{k y ) = v^—ky) are satisfied, it is possible to verify that rf(k y ) = 'y(-ky). This means 
the generation of Majorana bound state at the edge for k 2 > \k y \ > hi. A similar Majorana 
bound state also appears for A s > A t and /c 2 > \k y \ > k\. 

Unlike Majorana Fermions studied before, the present single branch of Majorana bound 
state is realized with time reversal symmetry. The TRI Majorana edge mode has the follow- 
ing three characteristics, a) It has a unique flat dispersion: To be consistent with the time- 
reversal invariance, the single branch of zero mode should be symmetric under k y — > —k y . 
Therefore, by taking into account the particle-hole symmetry as well, the flat dispersion is 
required. On the other hand, the conventional time-reversal breaking Majorana has a linear 
dispersion, b) The spin-orbit coupling is necessary to obtain the TRI Majorana edge mode. 
Without spin-orbit coupling, the TRI Majorana edge mode vanishes, c) The TRI Majorana 
edge mode is topologically stable under small deformations of the Hamiltonian (49). The 
topological stability is ensured by the topological invariant W{k y ). 

2. Z2 topological number and anisotropic response to the Zeeman magnetic field 

In this section, we would like to point out that the time-reversal invariant Majorana 
fermion is also characterized by another topological number taking a Z2 value. The Z2 
topological number explains the anisotropic response of the Majorana fermion to the Zeeman 
magnetic field found in Ref. 107 ' 108 . 

For a d xy + p-wave superconductor, the BdG Hamiltonian has the following symmetry, 



Regarding k y as a parameter, we can consider this as the particle-hole symmetry in one 
dimension. Thus the Z2 topological number can be introduced in a similar manner as 
shown in Ref. 67 . As seen from Appendix C, the Z2 topological number is given by (— l) u ^ 




■i-i 




(75) 
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TABLE IV: The zero energy ABSs of the semi-infinite d xy + p-wave Rashba superconductor on 
x > 0. In (a), we consider the spin-singlet dominant Cooper pair, A s > At, and in (b) the 
spin-triplet dominant one, At > A s . In the second and third columns, the chirality of each wave 
function ipf in (61) is summarized. As is explained in the text, the numbers of the zero 

energy ABSs with T = ±1 are determined from the chirality. For comparison, we also show the 
topological number W(k y ) given in (59) and (60). In both of the case (a) and (b), the results agree 
with the index theorem (17) excellently. 

with 

1 f°° 

u(k y ) = - / dk x A x (k). (76) 

JO 
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Here A x (k) is the "gauge field" defined by the bulk wave function \u n (k)}, 

A x (k) = tJ2Mk)\d kx u n (k)). (77) 

n 

Then the integral can be evaluated as 

(-1)"<*«> = sgn [(k 2 y /2m - - (\k y f] . (78) 

(See Eq.(C14).) From this, it is found that the Z 2 topological number is non-trivial, i.e. 
^—\y( k y) — —\ j n the region k\ < \k y \ < k 2 where the dispersionless Majorana fermion 
exists. Therefore, in addition to W(k y ), the Z 2 topological number (— l) 1/ ( fe y) also ensures 
the topological stability of the dispersionless Majorana fermion. 

The merit of the Z 2 topological invariant is evident if we apply a Zeeman magnetic 
field. In the presence of Zeeman magnetite field, the chiral symmetry (7) is broken since 
it is a combination of the particle-hole symmetry and the time-reversal symmetry, and the 
Zeeman magnetic field breaks the time-reversal invariance. This implies that the topological 
protection of the gapless state discussed in Sec. Ill does not work in the presence of a Zeeman 
magnetic field. In addition, the topological number W(k y ) cannot be defined without the 
chiral symmetry. 

On the other hand, the symmetry (75) survives even in the presence of the Zeeman 
magnetic field if we apply it in the ^/-direction. Indeed, for the d xy + p-wave superconductor, 
we can show that the BdG Hamiltonian 

y At(fe) -e(fc) + g(k) ■ a* + ^ B H y a* y ) 

in the presence of the Zeeman magnetic field H y in the y-direction satisfies (75). In a similar 
manner, we can introduce the Z 2 invariant, and we obtain 

(-1)^) = sgn [(k 2 y /2m - /i) 2 - (\k y ) 2 - ^ B H y f] . (80) 

(See.Eq. (C15).) Thus, the Z 2 number remains non-trivial in the region where the dis- 
persionless Majorana fermion exists. In fact, the region of k y in which the Z 2 topological 
number is non-trivial is extended in the presence of H y , which suggests that the magnetic 
field in the ^/-direction stabilizes the dispersionless Majorana fermion. This is a peculiar 
property that is not seen in other dispersionless ABSs. From the bulk-edge correspondence, 
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we can conclude that the time-reversal invariant dispersionless Majorana fermion survives 
even in the presence of the Zeeman magnetic field if we apply it in the y-direction. 

Here we notice that the Z 2 number is very sensitive to the direction of the Zeeman 
magnetic field we considered. While it is well defined in the presence of a magnetic field in 
the ^/-direction, it becomes meaningless if we apply a magnetic field in the other directions 
since in the latter case the symmetry (75) is broken. Therefore, the dispersionless Majorana 
fermion is also very sensitive to the direction of a magnetic field, which is consistent with 
the surface density of states calculated in Ref. 107 ' 108 . 

VI. INDEX THEOREMS: A PROOF OF BULK-EDGE CORRESPONDENCE 

Finally, we would like to prove the index theorems, (17), (18), (28) and (29). 

A. Strategy 

Before going into the details, we would like to outline our strategy. In order to prove the 
index theorems, we introduce an adiabatic parameter a in the Plank constant h, 

h = ah , (81) 

where h denotes an original value of the Plank constant. When a — > 0, we have a classical 
limit of h — > 0, and when a — 1, the system returns to the original. 

We will first consider the semi-classical limit S « 1 (a < 1). Using the WKB approxi- 
mation, we prove the index theorems in the following sections. ((94) and (96) in Sec. VI B, 
and (122) and (125) in Sec.VIC, respectively.) 

Then we adiabatically increase a — > 1 until h goes back to the original value Ho- From 
the argument in Sec. Ill, we notice here that Nq + ^ — cannot change in this process. 
Indeed, in order to change Nq + ^ — Nq~~^ , we need a continuum mode which closes the gap 
at the corresponding k\\, but we have a gap in the bulk from the assumption. A new state 
might appear near the boundary, but it should have a discrete spectrum for a fixed k\\ since 
it is localized near the boundary. Thus, the boundary state cannot change the value of 
Nq + ^ — as well. Therefore, the index theorems in the semi-classical limit remain to hold 
even when h goes back to the original value. 
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We would like to emphasis here that our strategy adapted here provides a general frame- 
work to prove the bulk-edge correspondence: For any topological state (of non-interacting 
systems), a mismatch of the topological number on the boundary results in a gap closing 
point near the boundary, in the classical limit (a = 0). Then, by using the WKB quanti- 
zation of the gap closing point, an gapless edge state is obtained semi-classically (a <C 1). 
From the existence of a gap in the bulk, the edge state is stable against an adiabatic change 
of a, thus we can conclude the existence of the edge state in the original theory (a = 1). 



B. 2x2 BdG Hamiltonian 

In this subsection, we prove (28) and (29). 

To make a boundary of a superconductor, let us introduce a confining potential V(x) 
illustrated in Fig. 8, 

^ A(fe) -e(k)-V(x) J 

We assume that the confining potential V(x) is steep enough near the edge. Now we will 
prove the index theorem (28) for the BdG Hamiltonian (82). 

First, we use the WKB approximation to count the zero energy ABSs. For this pur- 
pose, it is convenient to solve T-i2x2(k,x) 2 v = E 2 v instead of the original BdG equation 
%2x2(&, x)u = Eu. The zero energy states for these two equations are the same as is shown 
in Appendix. A. 

In the presence of the confining potential, T-i 2 x2(k 1 x) 2 is given by 

v (h ,2 ( (e(k) + V(x)) 2 + A(k) [V(x),A(k)} \ 

^2x2{k,x) = \ \ . (83) 

^ -[V(x),A(k)] (s(k) + V(x)) 2 + A(k) 2 J 

Here k x should be treated as —ihd x in the above. Then noting that 

r)V r)A r)V r)A 

[V(x), A(Jfe)] = ^[x, h x \° w + 0(h 2 ) = ih°^ a w + 0(h 2 ). (84) 



we have 



8V 8 A 

n 2x2 (k, x) 2 = [(e(fc) + V(x)) 2 + A(fc) 2 ]l 2x2 - hj^^a y + 0(h 2 ). (85) 
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(a) 



V(x) 




xO(~0) 



X 



(b) 




x<0 



x=x°(~0) x>0 



FIG. 8: (color online) (a) A semi-infinite superconductor on x > 0. The thick curve denotes 
the confining potential V(x), and the thin curve a gap of the system in the classical limit. If 
Eq. (87) is satisfied, the gap AE in the classical limit closes near the edge at x = 0. (b) The 
corresponding normal dispersion of electrons. Inside the superconductor (x > 0), the system 
supports the Fermi surface, but outside the superconductor (x < 0), the system becomes a band 
insulator. Correspondingly, AE is a superconducting gap for x > 0, and AE is a band gap (~ V{x)) 
for x < 0. 

In the classical limit (h — > 0) of the WKB approximation, the energy spectrum is given 
by the first term of the right hand side of (85), 



where k and x should be considered as c-numbers. For a fixed fen, we have a zero energy 
state if (k x ,x) = (k®,x°) satisfies 



Let us now take into account the leading order correction of h. Near the zero (h®,x°) 



E 2 = (e(k) + V(x)) 2 + A(k) 2 , 



(86) 




(87) 
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satisfying (87), we obtain 

e(fe) + V{x) 



dV(x) 



dx 



(x - x u ) + 



A(k) 



fdA(k)\ 



(k x - kl) + 



V dk x J kx=k o 

thus the first term of (85) is evaluated as a harmonic oscillator, 
Mk) + V(x)f + A(kf]= (|Q 2 ^z- 



2 

\ OK xy kx=k o 

From the Bohr-Sommerfeld quantization condition, it leads to 



(*. - KY + 



{e(k) + V(x)) 2 + A(k) 2 } = 2h 



dV 




dA 




dx 


x=x° 


dk x 


U —L.Q \ 
t^x—t^x 



n+-)+0(h 2 ), 



(88) 



(89) 



(90) 



with n — 0, 1, 2, • • • . The second term in (85) is evaluated as the expectation value for the 
WKB wave function, 



-h 



dV dA\ 



dxdk x / ° y V dx J x=x o \dk x J t 



fdA\ 



(91) 



b —bO 



n 2X 2(k,x) 2 = 2h 



dV 




OA 




dx 


x=x° 


dk x 


k —L.0 \ 



n+ D 12x2 - n (^ 



where we have used the fact that the WKB wave function has a sharp peak at x = x° and 
k x = k x . Combining (90) and (91), we obtain 

9 y) °v (92) 

For a semi-infinite superconductor on x > 0, the confining potential satisfies (dV/dx) x=x o < 
0. Therefore Eq.(92) yields a zero energy solution with chirality 7(= an eigenvalue of a y ) 
given by 

fdA 



7 = -sgn 



(93) 



Counting the zero energy states from all (k x ,x°)s satisfying (87), we obtain 



-J^sgn 



OA 



(94) 



,(+) «H 

W kx=k 

where the summation is taken for (k x ,x°) satisfying (87). 

Here we can show that (94) is nothing but the index theorem (28): First, we note that if 
and only if e(k) < 0, there exists a x satisfying e(k) + V(x) = 0. Thus we can rewrite (94) 
as 

fdA 



n, 



(+) 



i i n -■- - ^2 sgn 

A(fcg,fc,|)=0,e(fcg,fe||)<0 



\9k x J kx=k o ' 



(95) 
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As is shown in Sec.B 1, the right hand side of the above equation is equal to —w(k\\), thus 
we have the index theorem 

nJ -) -n { +) = w(k\\). (96) 

in the leading order in the WKB approximation. 

While we have derived Eq. (96) in the leading order of the WKB approximation, we can 
find that no further correction exists: As was discussed in Sec. Ill, any small perturbation 
preserving the chiral symmetry cannot change — unless the bulk energy gap closes. 
Therefore, a higher order correction might change each of or Uq \ but it cannot change 
their difference — \ This means that Eq.(96) is exact. (See also the discussions in 
Sec. VI A.) 

In a similar manner, we can prove the index theorem (29) for a semi-infinite supercon- 
ductor on x < 0. Following the same argument above, we obtain (92) again. However, 
in comparison with the semi-infinite superconductor on x > 0, the sign of (dV/dx) x=x o is 
reversed for a semi-infinite superconductor on x < 0, so the chirality 7 of the zero energy 
state is also reversed 

7 = sgn (w) • (97) 



Thus we obtain 



(+) „(-) ^ 1 dA 



^{bk)^' (98) 

where the summation is taken for (k x ,x°) satisfying (87). This leads to 

n ( +) -n^ =w(k ll ), (99) 
which is exact again for the same reason above. 



C. General case 



Now, we prove the index theorems (17) and (18) for general time-reversal invariant BdG 
Hamiltonian. 

First let us consider the index theorem (17) for a semi-infinite superconductor on x > 0. 
To realize the semi-infinite superconductor, we introduce a confining potential V(x) in a 
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similar manner as Sec. VI B (See also Fig.8), 

H(k)^nk,x)=( i(kU ' + V(X)!i *" A<feW |, (100) 

y At(fe)„„, -E T (-k) ac ,-V(x)S ac 

where V(x) = inside the superconductor (x > 0). We suppose a sharp edge where the 
confining potential V(x) is very steep near the edge (x — 0). 

In the classical limit h — > 0, k x and x commute with each other, thus we can treat them 
as c-numbers. Then, considering i as a parameter, we define the winding number W(k\\,x) 
in a manner similar to Eq.(ll), 

W{k h x) = — — / dk ± tr \m- l (k,x)d k± H(k,x)} 
Arci J 

= — Im / dk ± d k± In det q(k, x) , (101) 

27T [J 

where q{k,x) = i[S(k) + V{x)]U^ — A(fe) with U in Eq.(4), and = k x . The line integral 
in the above is defined in the same way as Eq.(ll). 

Here we find that W(k\\,x) is identical to W(k\\) inside the superconductor x > 0, 

W(k h x) = W(k\\), for x > 0, (x ^ 0), (102) 

since V(x) = there. We also find that W(k±,x) becomes zero outside the superconductor, 

W(fe||,a;) = for x < 0, (x ^ 0), (103) 

since we have a vacuum there. Therefore, if the bulk superconductor supports a non-zero 
winding number W(k\\) ^ 0, then the value of W(k\\,x) must be changed near the edge 
(x — 0). This immediately implies that when W(fe|i) ^ 0, a gap of the system closes near 
the edge in the classical limit: Since W(k\\,x) is a topological number, it changes only when 
the line integral is ill-defined. Then this occurs only when det q(k, x) = 0, which means a 
gap closing of the system in the classical limit. 

To examine how W{k\\,x) changes, we diagonalize q(k,x) as 

q(k,x) = A(k,x)A(k,x)B' ( (k,x), 
q\k, x) = B(k, x)A*(k, x)A\k, x), (104) 

where A(k,x) = diag(Ai, A2, • • • , Ajv) with Aj the eigenvalue of q(k,x), and A(k,x) and 
B{k,x) are N x N matrices satisfying 

A(k, x)B^(k, x) = l NxN . (105) 
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(Here we assume that q(k,x) is an N x N matrix. See Appendix E for details.) Since 
Eq.(104) yields that detq(k,x) = detA(k,x) = Y[ n A a (k,x), the winding number W{k,x) 
is recast into 



w ( k h x ) = 7^: Im 



N 
n=l J 



dkj_dk ± In X n (k, x) 



(106) 



Thus the line integral of W(k\\,x) is ill-defined, when some of the eigenvalues \ n (k,x) have 
a zero. 

Now suppose that \ n (k,x) has a zero at (x, k x ) = (x°, k x ), 



ae 



03 



{x-x°)+ 1 e i5 {k x -k° x ) + 



= ae^ [(x - x°) + (^/a)e t{5 -^(k x - k° x )] + 



(107) 



where a, (3, 7 and 5 are real constants. We notice here that a is very large since it is 
proportional to (dV/dx) x=x o and V(x) is steep near the edge. Thus we can neglect the real 
part of ('j/ay^ s ~ l3 ' ) (k x — and obtain 

A n = ae lfi [{x - x°) + ir]{k x - k° x )] +■■■ , (108) 

with r] = (7/a) sin(5 — (5). This zero changes the value of the winding number W{k\\,x) as 
follows. Consider the line integrals near the zero along x = x % illustrated in Fig 9, 

/ dk x d k Jn\ n (k,x) (2 = 1,2) (109) 

J x=x l 

with x 1 < x° < x 2 . The difference of these line integrals is evaluated as 



dk x d kx lnA n (fe,x) - / dk x d kx In A„(fe, x) 

= j> dk x d kx \n\ n (k,x) 
= 27nsgn[?7], 



(110) 



thus from Eq.(106), we have 



W{k h x 2 ) - W(k h x l ) = sga[rj\. 



(Ill) 



In a similar manner, we can consider other zeros of the eigenvalues of q(k,x). Then 
summing up all contribution of these zeros k x ^) of A„(j) (i — 1, • • • , M), we have 

M 

H/(fc||) = ^sgn[^], (112) 



i=l 
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FIG. 9: (color online) A zero (x°,k®) of \ n near the edge of a semi-infinite superconductor on 
x > 0. 



where rji is given by 

A n(i) = [(x - x ^) + irji(k x - + • ■ ■ , (113) 

with real functions cti and fa. Since detT-L(k,x) around the zero fcjj!^) is evaluated as 

detH(k,x) oc det g(fc, x)q*(k, x) 

= n|An(fc,x)| 2 

~ [(x-x^) 2 + i^(A«-A?fl) a ], (114) 

the classical spectrum around each zero is a harmonic oscillator. 

Let us now take into account the quantum corrections. In the leading order correction 
of h, the classical spectrum of a harmonic oscillator around each zero is quantized. To see 
whether the zero energy state survives or not after the quantization, let us examine the BdG 
equation for the zero energy state, 

H(k 1 x)u = 0i (115) 

which is equivalent to 

q{k,x)<p = 0, g t (fc,x)^o = 0, (116) 

with 




(117) 



Substituting the following expression 



0o = a n{i) (k U,k h x°V)f(x), ij = b n(l) {kf\k h x^)g{x), 



(118) 



into Eq.(116), we find that 



[a; - x 0{i) + irji(k x - fc£ (i) )] f(x) = 0, 
[ x - x °W-ir h (k x -kW)]g(x)=0, 



(119) 



around the zero (x°^\k x ^). Here a n (k,x) and b n (k,x) are the eigenfunctions of q(k,x) 
and q'(k,x) in the classical limit, defined in Appendix E. Replacing k x with —ihd x , we can 
easily solve Eq.(119), which reads 



f(x) = /oexp 
g(x) = g exp 



' k 0(i) 

I— — X 



h 



cxp 



' k m 

1——X 

n 



cxp 



2rjih 

(x-x ^) 2 
2i]ih 



(120) 



where fo and go are constants. Since the normalization condition requires that fo is zero if 
r/i < 0, and g is zero if rji > 0, the zero energy state has r = 1 for rji < and T = — 1 for 
rji > 0, respectively (We have used Eq.(8) when applying T to u in Eq.(117).) In other 
words, in the leading order of the WKB quantization, there exists a zero energy state for 
each zero (x°( l \ k x ^), and its chirality T is given by 



T = -sgn^]. 



(121) 



From (112), this leads to 



W(kn) = N ( 



(-) 



N, 



(+) 



122) 



From the same argument as in Sec. VI B, we can conclude that no higher order correction of 
h exists in Eq.(122) again. Then using the argument in Sec. VI A, we conclude that Eq.(122) 
is an exact result. 

In a similar manner, we can derive (18) for a semi-infinite superconductor on x < 0. As 
illustrated in Fig. 10, for a zero (x°, k®) of \ n near the edge x = 0, we have again Eq.(llO), 



L 



dk x d kx In A n 



dk x d kx lnA n = 27risgn[r?], 



(123) 
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FIG. 10: (color online) A zero (x°,k®) of A n near the edge of a semi-infinite superconductor on 
x < 0. 

with a; 1 < a; < x 2 , but in contrast to the previous case, the line integral along x = x 2 is 
located at an outer region of the superconductor. Thus the sign of the right hand side of 
Eq.(112) is reversed, 

W(k {] ) = -J2sgn[ Vi ], (124) 

i 

and the following index theorem is obtained 

W(k {] ) =iV (+) -iV H . (125) 
From the argument in Sec. VI A, this result holds again beyond the WKB approximation. 

VII. CONCLUSION 

In this paper, we have discussed the topology of the ABS with flat dispersion at zero 
energy which appears on a boundary of time-reversal invariant superconductors. Using the 
symmetry of the BdG Hamiltonian, we have introduced the topological numbers W(fc||) and 
w(k»), and from the bulk-edge correspondence, topological criteria for the dispersionless 
ABS were obtained, which are summarized as the index theorems (17), (18), (28) and (29). 
We have shown that the index theorems correctly predict the existence of the dispersionless 
ABSs. It has been also clarified that sign change of the gap function is directly related to 
our topological criterion for superconductors preserving S z with simple Fermi surfaces. As 
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concrete examples, we have also discussed (i) d^-wave superconductor, (ii) p^-wave super- 
conductor, and (iii) d xy +p-wave superconductor. All the examples confirms our topological 
criteria excellently. In the last part of this paper, we provide a general framework to certify 
the bulk-edge correspondence, and we prove the index theorems. 

While we have considered only superconductors in this paper, the index theorems we 
found can apply to other systems. For example, graphene in the nearest neighbor tight- 
binding model has chiral symmetry (or sublattice symmetry) similar to (23) 130 , thus the 
index theorems (28) and (29) hold. In general, our index theorems apply to any non- 
interacting Hamiltonian with chiral symmetry. 
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Appendix A: basic properties of system with chiral symmetry 

Here we summarize the properties of the system with chiral symmetry. Suppose that the 
Hamiltonian % has the chiral symmetry Y 

{n,r} = o, r 2 = i. (Ai) 

Then consider the BdG equation 

U\u E ) = E\u E ) (A2) 

First, we would like to show that there exists a one to one correspondence between the 
solution of (A2) and that of the following equation, 

n 2 \v E 2) = E 2 \v E 2). (A3) 

Since it is evident that the solution \u E ) of (A2) satisfies (A3), here we show only the 
reverse. Namely, we can construct the solution \u E ) of (A2) from the solution \v E t) of (A3): 
For E 2 ^ 0, if (H + E)\v E 2) ^ 0, the solution \u E ) of (A2) is obtained by 

\u E ) = c{U + E)\v m ) (A4) 
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with a normalization constant c. Then if ("H + E)\ve*) = 0, the solution is given by 

\u E )=T\v E 2). (A5) 

It is easily found that \ue) in (A4) and (A5) has a nonzero norm and satisfies (A2). For 
E 2 = 0, we have 

n 2 \v )=0. (A6) 

Since this equation implies that 

( Vo \n 2 \v ) = \\n\v )\\ 2 = o, (A7) 

we obtain 

H\v ) = 0. (A8) 

Therefore, the zero energy solution \u ) is obtained by \u ) = \v ). Since there exists a 
one to one correspondence between the solution of (A2) and that of (A3), we examine (A3) 
instead of (A2) in the following. 

Because V, 2 and T commute with each other, we can take the solution of (A3) as the 
eigen state of T, simultaneously. Then, it is found that for E 2 ^ 0, the solution |t>+ 2 ) 
satisfying r|t>+ 2 ) = \v~^ 2 ) is constructed from the solution \v^ 2 ) satisfying r|i>~ 2 ) = — \ v e 2 ) 
by multiplying T-L from the left, 

\v+ 2 ) = c"H\v- 2 ), (A9) 

where d is a normalization constant. In a similar manner, we can construct \v~^ 2 ) from \Vg 2 )- 
Therefore, for E 2 7^ 0, the solution of (A3) with the chirality T = 1 is always paired with 
the solution with T = — 1. 

On the other hand, for E 2 = 0, the solution does not form a pair in general. Indeed, as 
was shown in the above, the solution \v ) with E 2 = satisfies V,\vq) = 0, thus Eq.(A9) 
leads to \vq) = 0. Thus we do not obtain the paired state in this case. 
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Appendix B: Derivation of useful formulas 



1. Formulas for w(k y ) 



We first evaluate the integral (26) 



w(k v ) = — Im 



J dk x d k Jnq(k) , (Bl) 
with q(k) = —ie(k) — A(k). To calculate this integral, it is convenient to rewrite this as 

w(k y ) = ~^J dk x e ab m a (k)d k m b (k), (B2) 

where 

e(fc) ... A(fc) 



mi(fc) 



m 2 (k) 



(B3) 



y/e 2 (k) + A 2 (fc)' y/e 2 (k) + A 2 (k)' 

Then we apply a technique used in Ref. 123 ' 126 . Since the integral (B2) is a topological 
number, it can not change its value even if we rescale A(k) as aA(k) (a < 1), 

e(Jfe) /n aA(fe) 



mi(fe) — >■ 



m 2 (k) 



(B4) 



V£ 2 (fc) + a 2 A 2 (fe) ^(fc) + a 2 A 2 (fc) 

Then when o<l, except the neighborhoods of the zero of m a (k) becomes a constant 



mi(fc) ~ ±1, m 2 (k) ~ 0. 



(B5) 



Thus only the neighborhood of the zero of e(k) contributes the integral. Expanding e(k) 
and A(k) around the satisfying e(k®, k y ) = 0, 



e(fc) = d kx e(k° x , k y )(k x -k° x ) + ---, A(k) = A(k° x , k y ) + 



(B6) 



we estimate the contribution near the as 



-sgn[d kx e(k° x , k y )}sgnA(k° x , ky). 
Summing up the contribution from fc°'s, we obtain Eq.(27) 



(B7) 



(B8) 



e(kO, ky )=Q 
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FIG. 11: Examples of A(fc) as a function of k x . The symbols db represent the sign of d kx A(k x ) at 
k x = k x with A(k x ) = 0. (a) A(fe) in a lattice model. From the periodicity in k x , Eq.(BlO) holds, 
(b) and (c) A(fc) in a continuum model. Here we regulate A(fc) as A(fc) — > at /c x = ±oo. In the 
case (c), we need to take into account the contribution at k x = oo to obtain Eq.(BlO). 



In a similar manner, we also obtain the following formula by exchanging e(k) with A(fc) 
in the above, 

w ( k v) = ~\ E s § n [VWi *v)] s S n ^°> K)- (B9) 

A(fcg,fc y )=0 

As illustrated in Fig. 11, the following identity holds in general, 

E sgn^A^A;,)] =0. (BIO) 

A(feO,fe !/ )=0 

For a lattice model, this equation comes from the periodicity in k x , and for a continuum 
model, this holds from the regularization in which A(fc) at fc^ = oo is identified with that 
at k x = — oo. Since the left hand side of Eq.(BlO) is rewritten as 

E sgn[^A(A£A)] 

A(fcg,fcy)=0 

E sgn[d kx A(klk y )] 

A(kO,k y )=0,e(kO,ky)>0 

+ E sgn[d kx A(k° x ,k y )], (Bll) 

A(fcg,k,,)=0,6(ft2,*»)<0 

Eq.(BlO) leads to 

E sgn[d kx A(k x ,k y )}=- E sgn[d kx A(k x ,k y )} (B12) 

A(fcg,fc !/ )=O, £ (A ; g,fej / )>0 A(fcg,fc J ,)=0, e (feg,fc y )<Q 

42 



From this equation, Eq.(B9) is recast into 

1 



w(k y ) = -- sgn[d kx A(k° x ,k y )} 

A(k x ,k y )=0,e(k x ,k y )>0 

sgn[d kx A(k° x ,k y )} 

A(k x ,k y )=0,e(k x ,k y )<0 

sgn[d kx A(k° x ,k y )]. 



1 

+ 2 



A(k x! k y )=0,sgne(k x ,k y )<0 

We will use this formula in Sec. VI B. 



(B13) 



2. Formulas for W(k y ) 



Here we evaluate the integral (51) 



W(k y ) = -Im 



dk x dk x lndetg(fc) 



with 



q(k) = [e(k) — iip(k)]a y + [g(k) — id(k)} ■ cra y 



Since the determinant of q(k) is given by 



(B14) 



(B15) 



detq(k) = -e 2 {k) + g 2 (k) + ^ 2 (fc) - d 2 (k) - 2i [g(k) ■ d(k) - e(k)^(k)] , (B16) 

the above integral is recast into the same form as (B2) with 

2(g(k) ■ d(k) - e(k)^(k)) 



mi(fc) = 
m 2 (k) = 



y/W(k) - g 2 (k) - ^ 2 (fc) + d 2 (k)] 2 + A[g{k) • djk) - e(k)^(k)] 2 ' 
e 2 {k) - g 2 {k) - ij 2 {k) + d 2 {k) 



(BIT) 



^[s 2 (k) - g 2 (k) - ^ 2 {k) + d 2 (k)] 2 + 4[g(k) ■ d(k) - e(k)^(k)] 2 ' 

For a weak pairing superconductor, the energy scale of the gap functions ip(k) and \d(k)\ 
is much smaller than that of e(k) and \g(k)\. Thus we can rescale ip(k) and d(k) as a^p(k) 
and ad(k) with a positive small constant a without changing the value of W(k y ). Then for 
a <C 1, it is found that only the neighborhood of the zeros of e 2 (k) — g 2 (k) contributes to 
W(k y ). As a result, in a similar manner to Appendix B 1, we obtain the formula (52). 
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Appendix C: Zi topological number 



First consider the quasi-particle wave function \u n (k)) with 

n(k)\u n (k x , k y )) = E n (k)\u n (k x , k y )). 



(CI) 



From the symmetry (75), we can say that if \u n (k)) is a quasiparticle state with posi- 
tive energy E n {k) > 0, then C'\u* n (—k x ,k y )) is a quasiparticle state with negative energy 
—E n (—k x , k y ) < 0. (See Eq.(75) for the definition of C.) Thus we use a positive (negative) 
n for \u n (k x , k y )) to represent a positive (negative) energy quasiparticle state, and set 



\ u -n(k x , k y )) — C \u n ( k x ,k y )). 

For the ground state in a superconductor, the negative states are occupied. 
Now we introduce the following "gauge field" A x ± \k), 



(C2) 



(C3) 



A x ± \k)= i J2Mk)\d k u n (k)}. 
From (C2), the gauge field A x ± \k) satisfies 

A x + \k X ,ky)=A x ~\~k x ,ky). (C4) 

We also find that A x (k) = A x +) (k) + A { x \k) is given by a total derivative of a function. To 
see this, we rewrite \u n (k)) in a component, 

f<(k)\ 



\u n {k)) 



<{k) 
<{k) 

and introduce the following 4x4 unitary matrix V(k) 

V mn {k) =<(fc), (m= 1,2,3,4, n = 2,l,-l,-2). 
Then A x (k) is written as 

A x (k) = iti{V\k)d kx V{k)) = id kx lndety(fc). 



(C5) 



(C6) 



(C7) 



In a similar manner as in Ref. 67 , we introduce the Z 2 topological number (— l)^) as 



^(J^y) — / dk x A x ^{k x ,ky 



(C8) 
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Here we regulate the gap function as A(k) — > far apart from the Fermi surface to define 
the topological number. Eqs. (C4) and (C7) lead to 

If 1 f°° 

^(ky) — / dk x A x ^(k x ,k y ) H / dk x A x ^(k x ,k y ) 

K J-oo 7T J 

1 f 00 

— — dk x \_A X + ^ (k x , k y ) + A x \k x , k y )\ 
71 Jo 

(C9) 



7T 



dety(oo, fcy) 
detF (0, fcj,) 



thus, we obtain 



(-!)"<*«) = detV( °'^. (CIO) 
1 ; dety(oo,y 1 ; 



We can see that the right hand side of (CIO) indeed takes a Z 2 value: From Eq.(C2), V(0, k y ) 
is given by 

V(0,k y ) = (\u 2 (0,k y )),\ Ul (0,k y )),Cyi(0,k y )),C'^*(0,k y ))), (Cll) 
thus it satisfies 

C'V(0,k y ) = (Cy 2 (0,k y )),C'\ Ul (0,k y )),\ul(0,k y )),\u* 2 (0,k y ))). (C12) 

Calculating the determinants of both sides, we obtain detV(0, fc y ) = detV*(0,k y ). This 
implies detl / (0, fcj,) = ±1 since V(0, k y ) is a unitary matrix. In a similar manner, we obtain 
dety(oo, k y ) = ±1 from the regularization in which k x = oo and fc^ = — oo are identified. 
So the right hand side of (CIO) takes ±1. 

Using the explicit form of the quasiparticle wave functions of (49), we obtain 

detV(0, ky) = -sgn[(£; 2 /2m - /i) 2 - {Xk y ) 2 ] (C13) 

for the d xy + p-wave superconductor. Furthermore, due to the regularization of A(k), we 
find that H(k) is almost diagonal at k x = oo. Thus we have |u2(oo, k y )) ~ T (l, 0, 0, 0), and 
|tii(oo, fcj,)) ~ T (0, 1,0,0), which implies detV r (oo, k y ) = — 1. Therefore, the Z 2 topological 
number is evaluated as 

{-iy^ = sgn[(A; 2 /2m - /i) 2 - (Xky) 2 }. (C14) 

Note that there are two differences between the above Z 2 topological number and that in 
Ref. 67 . First, the Z 2 topological number (C8) is defined for any k y , while the Z 2 topological 
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number in Ref. is defined only for the time- reversal invariant k y . This difference comes 
from the difference in the symmetry we used. In contrast to the Z 2 topological number in 
Ref. 67 , where the particle-hole symmetry (6) was used, in order to define (C8), we assume 
the one- dimensional particle-hole symmetry (75) which requires a special symmetry of gap 
function such as d xy +p-wave pairings. Since the Z 2 topological number in Ref. 67 is available 
only for the discrete values of k y , it cannot explain the existence of the ABS with flat 
dispersion, although the Z 2 topological number (C8) can. Second, the line integral in (C8) 
performed from k x = to k y = oo as we consider the continuum model of the d xy + p-wave 
superconductor. 

Even in the presence of the Zeeman magnetic field H y in the y-direction, the symmetry 
(75) persists in the d xy + p-wave superconductor. Thus, the Z 2 topological number can be 
defined in the same manner as above. The resultant Z 2 number is 

(-1)^) = sgn[(A; 2 /2m - /i) 2 - (AA^) 2 - {^H y f\. (CIS) 

Appendix D: Relation to the odd- frequency pairing amplitude generated at the 
surface 

In this appendix, we would like to discuss the relevance of Andreev bound state and odd- 
frequency pairing function in more detail 49 ' 51 ' 52 . It has been clarified that inhomogeneity of 
superconductivity can induce odd-frequency pairing function 49 ' 51 ' 52 ' 131-134 . For this purpose, 
it is better to use quasi classical Green's function, where atomic scale oscillation is removed 
out 135-143 . Here, we assume Cooper pair with S z — for simplicity. 

First, we consider two-dimensional semi-infinite superconductor in x > where flat sur- 
face is located at x — 0. In the following discussion, to express the symmetry of the frequency 
dependence of Cooper pair explicitly, we use Matsubara frequency. Quasiclassical Green's 
function can be written as 

9+ = _ 
9- = _ 

V- 



u 

-9+ 



1 + D + F: 



2iF. 



+ 



2iD + -(1 + D+F+) 



/{1-D+F+), (Dl) 



f- 



1 + D_F_ 
-2iF_ 



1 + 



/(1-£>_F_), (D2) 
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by using Riccati parameterization where g and /(/) are normal and anomalous Green's 
functions, respectively. The subscript +(— ) denotes the right and left going quasiparticles. 
The above Riccati parameters D± and F± obey the following equations, 

v Fx d x D± = -A±(l - D 2 ± ) + 2u n D ± , 

v Fx d x F ± = -A±(l - Fl) - 2u n F±. (D3) 
The bulk solution of D± for x > can be written as 

D± = A± (D4) 

Here, A± is a pair potential felt by quasiparticle for each trajectory. v Fx denotes the x- 
component of Fermi velocity with v Fx > 0. For bulk d^-wave superconductor, it is given 
by 

A ± = ^A , 

for k x > 0. On the other hand, the corresponding quantity for p^-wave, it is given by 

* i * 

"=A , 



Vfc 2 

for k x > 0. 

At the surface x = 0, F + = —D_ and F_ = — D + are satisfied due to the specular 
reflection. Also in the presence of ABS, D + = —D_ are satisfied. Thus, we obtain 



'♦ = '- = i^r < D5 » 



Since D + (— ui n ) = l/D + (u) n ) is satisfied, we can show that /+ and /_ is an odd-function of 
Matsubara frequency to n . This means that a purely odd- frequency pairing state is realized 
at the surface in these two cases. 

In order to clarify the parity of the pairing, we define pairing amplitude which is available 
both for positive and negative k x . We write explicitly the k dependence of the pairing 
amplitude. We introduce f(k x , k y ) as f(k x , k y ) = f+(k x , k y ) and f(—k x , k y ) = f-(k x , k y ) for 
k x > 0. For even-parity pairing, f(k x ,k y ) = f(—k x ,—k y ) is satisfied. On the other hand, 
for odd-parity case, f(k x , k y ) = —f(—k x , —k y ) is satisfied. Here, let us discuss the parity of 
the odd-frequency pairings obtained above. For d xy -w&ve pair potential, since f-(k x ,k y ) = 
—f~(k x ,—k y ) is satisfied, it is straightforward to derive that f(k x ,k y ) = —f(—k x ,—k y ) 
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is satisfied. This means the realization of odd-parity pairing in this case. For p x -wave 
pair potential, since f-(k x ,k y ) = f-(k x , —k y ) is satisfied, it is easy to see that f(k x ,k y ) = 
f(—k x , —k y ) is satisfied. This means the generation of even-parity pairing in this case. 

To see concrete form of f(k x , k y ), we neglect the spatial dependence of A±. Then, f± is 
given by 

for d xy -wave pairing and 

/+ = /. = (D7) 
for p^-wave pairing. The resulting f(k x , k y ) is given by 

, \ _ ^sgn(/c :r )Ao k x k y _ isgn(k y )A \k x \\k y \ 
J\k x , ky) — - — - - - , (D8) 

for d xy -wave pair potential, and 

/(*,, ky) = ^ k ^ j^ = (D9) 

for p^-wave pairing. It is evident that odd-frequency odd-parity pairing in the former case 
and odd-frequency even-parity pairing in the latter case are realized. This is consistent 
with the fact that the spin of Cooper pairing is conserved. The difference of the parity 
results in a remarkable difference when we consider proximity effect into DN attached to 
superconductor 49 ' 51 ' 52 . 

In DN, only s-wave even parity pairing is possible due to impurity scattering. Thus, 
odd-parity odd-frequency pairing amplitude cannot penetrate into DN. Thus, for d xy -wave 
superconductor, ABS cannot enter into DN since it is expressed by odd-frequency spin- 
singlet odd-parity state. On the other hand, for p^-wave superconductor, ABS can enter 
into DN since it is expressed by odd-frequency spin-singlet even-parity state including s- 
wave channel. Thus, it has been shown that proximity effect of d xy -wme superconductor 
and p^-wave superconductor is completely different 44-48 . 

Next, we consider a two-dimensional semi-infinite superconductor in x < where flat 
surface is located at x — 0. Quasiclassical Green's function can be written as 




1 + D + F + -2iD + 
-2iF + -(1 + D + F + ) 



/(1-D + F + ), (D10) 
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, 9- f- \ / 1 + D-F- 2iF^ 
9- = = , . \ (1-D.F-), Dll 

/- -9- J \ 2tD_ -(1 + D_F_ 

by using Riccati parameterization. The subscript +(— ) again denotes the right and left 
going quasiparticles. Following the similar discussion in the case where two-dimensional 
semi-infinite superconductor is located on x > 0, we can use the relation D + = — F_ and 
D + = —F + . In the presence of Andreev bound state (ABS), D + = —D_ are satisfied. Then, 
we can write 



-2iD 

i-d; 

If we neglect the spatial dependence of A±, f± is given by 



/+ = i FiT - (D12) 



/+ = / _ = _^w»_ (D;I:5) 

Uj n ft/ 

for spin-singlet d xy -wave pairing and 



/+ = /- = — "7=, (D14) 

for spin-triplet p^-wave pairing. Following similar discussions below Eq. (D6), f(k x ,k y ) is 
given by 

f(l , x isgn(k x )A Q fca-fcy isgn(fc y )A \k x \\k y \ 
for spin-singlet d xy -wave and 

/(*.,*,) = -SS^ 4= = -^M (D16) 
for spin-triplet p^-wave. 

Appendix E: diagonalization of q(k, x) 

Let us consider an N x A matrix x) with eigenvalues Ai, • • • , Ajv- Then we have two 
kinds of eigenvectors : The first one is right eigenvectors, 

q(k, x)a n (k, x) = X n (k, x)a n (k, x), (El) 

and the second one is left eigenvectors 

bl(k,x)q(k,x) = bi(k,x)X n (k,x), (E2) 
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with the normalization condition 



bl(k, x)a m (k, x) = S mn . (E3) 

Note that a n (k,x) ^ b n (k,x) in general unless q(k,x) is a hermitian matrix. Eqs.(El) and 
(E2) lead to 

q(k,x)A(k,x) = A{k, x)A(k, x), 
q\k,x)B(k,x) = B(k,x)A*(k,x), (E4) 

with 

A(k,x) = (ai(fe, x), • • ■ ,a N (k,x)), 
B(k,x) = (6i(fe,a;), • • • ,b N (k,x)), 

A(k,x) = diag(Ai(fe,x), • • • , \ N (k,x)). (E5) 
Using the normalization condition (E3) which reads 

B f (fc, x)A(k, x) = A\k, x)B(k, x) = l NxN , (E6) 

we obtain 

q(k,x) = A(k,x)A(k,x)B\k,x), 
q\k, x) = B(k, x)A*(k, x)A\k, x). (E7) 
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